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Abstract Symbiotic stars are interacting binary systems with the longest orbital periods. They 
are typically formed by a white dwarf, a red giant and a nebula. These objects are natu¬ 
ral astrophysical laboratories for studying the evolution of binaries. Current estimates of the 
population of Milky Way symbiotic stars vary from 3000 up to 400000. However, the cur¬ 
rent census is less than 300. The Large sky Area Multi-Object fiber Spectroscopic Telescope 
(LAMOST) survey can obtain hundreds of thousands of stellar spectra per year, providing a 
good opportunity to search for new symbiotic stars. In this work we detect 4 of such binaries 
among 4,147,802 spectra released by the LAMOST, of which two are new identifications. 

The first is LAMOST J12280490-014825.7, considered to be an S-type halo symbiotic star. 

The second is LAMOST J202629.80+423652.0, a D-type symbiotic star. 

Key words: binaries: symbiotic-stars, optical spectra, LAMOST 

1 INTRODUCTION 

Symbiotic stars are interacting binary systems with the longest orbital periods and a multicomponent struc¬ 
ture, that is an evolved cool giant and an accreting hot, luminous companion (usually a white dwarf, WD; 
some might be a neutron star or a black hole) surrounded by a dense ionized symbiotic nebula. Therefore, 
symbiotic stars have composite spectra. 

Depending on the nature of the cool giant, symbiotic stars are generally classified into three types: 
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T e s ~ 3000 — 4000 K. They are the most numerous (~ 80%) and their orbital periods are of the order of a 
few years. D-types contain Mira-type variables and warm (T ~ 1000 K) dust shells. Their orbital periods 
are larger than 15 years (Mikolajewska 2003). D'-types contain F- or G-type cool giant stars surrounded 
by a dust shell with a temperature not higher than 500 K (Allen 1984b). There are only 8 D'-types out 
of 188 symbiotic stars in the Belczynski et al. (2000) catalogue. The presence of both an evolved giant 
that suffers from heavy mass loss and a hot component copious in ionizing photons results in rich and 
luminous circumstellar environment. In particular, very different surroundings are expected, including both 
neutral and ionized regions, dust forming regions, accretion disks, interacting winds and jets, which makes 
symbiotic stars a natural laboratory for studying the binary star interaction and evolution (Corradi et al. 
2003, and references therein). 

Symbiotic stars are proposed as probable progenitors of type la supernovae (SNe la) (Hachisu et al. 
1999; Chen et al. 201 1), as rapid mass accretion from the cool giant component may lead the WD to reach 
the Chandrasekhar mass limit (Han & Podsiadlowski 2004; Wang & Han 2012). SNe la have been suc¬ 
cessfully used as a cosmological distance indicator (Riess et al. 1998; Perlmutter et al. 1999), however, the 
exact nature of the progenitors of SNe la is still not clear. SNe la PIF 1 lkx, with a broad Ha line in its late 
spectrum, may result from a symbiotic star (Dilday et al. 2012). Large numbers of symbiotic stars can help 
us to study the progenitors of SNe la. Hence, searching for new symbiotic stars is of extreme importance. 

Current estimates for the total Galactic symbiotic population vary considerably from 3000 (Allen 
1984b), 30000 (Kenyon et al. 1993) to 3-4xl0 5 (Magrini et al. 2003). However, there are less than 300 
Galactic symbiotic stars known (Belczynski et al. 2000; Miszalski et al. 2013; Miszalski & Mikolajewska 
2014; Rodrfguez-Flores et al. 2014). This incompleteness is due to the lack of deep all-sky surveys at wave¬ 
lengths where symbiotic stars can be efficiently selected against other stellar sources (Belczynski et al. 
2000; Miszalski & Mikolajewska 2014; Mikolajewska et al. 2014). Most of the known symbiotic stars 
have been discovered by large Ha surveys, e.g. the INT Photometric Ha Survey and the AAO/UKST 
SuperCOSMOS Ha Survey, and then confirmed by the deep spectroscopy and long-term I-band light curves 
of the candidates. The Large sky Area Multi-Object fiber Spectroscopic Telescope (LAMOST) survey can 
obtain hundreds of thousands of stellar spectra per year, providing a good opportunity to find new symbiotic 
stars. Since no effort has been made from the large spectroscopic observation survey, we devise a selection 
method to find out symbiotic stars from the data released by LAMOST (including DR1, DR2 and DR3). As 
a result, we detect 4 symbiotic star spectra, of which two are new identifications, the other two have been 
previously identified by Downes & Keyes (1988). 

The structure of the paper is as follows. Section 2 briefly introduces the LAMOST survey and our 
method to select symbiotic stars. Section 3 describes the 4 symbiotic stars we identify. Section 4 is the 
discussion. We conclude in Section 5. 

2 SPECTROSCOPIC OBSERVATIONS AND TARGET SELECTION 
2.1 The LAMOST Survey 


LAMOST is a quasi-meridian reflecting Schmidt 4.0 m telescope located in Xinglong Station of National 
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Table 1 Data of LAMOST. 


Data 

Total 

Star 

Galaxy 

QSO 

Unknown 

DR1 

2,204,860 

1,944,406 

12,082 

5,017 

243,355 

DR2 

1,443,843 

1,240,792 

23,967 

4,251 

174,833 

DR3 

499,099 

437,299 

5,229 

574 

55,997 


Notes: The number of spectra, which are classified as stars, galaxies, QSOs and unknowns by the 
LAMOST pipeline. DR1 includes both 1,487,200 spectra obtained during the first year of the general sur¬ 
vey, and 717,660 spectra which were observed during the pilot survey. For more details see the website 
http://dr.lamost.org/ucenter/mytoken. 


scope is equipped with 16 low-resolution spectrographs, 32 CCDs and 4,000 fibers. Each spectrograph is 
fed with the light from 250 fibers and approximately covering a wavelength range 3700 — 9100 A at a 
resolving power of about 1800. With the capability of its large field of view and a large number of targets, 
LAMOST telescope dedicates for spectral survey of celestial objects in the whole northern sky (Cui et al. 
2012; Zhao et al. 2012). 

The current data of LAMOST contain data release 1 (DR1), data release 2 (DR2) and data release 
3 (DR3). The number of DR1, DR2 and DR3 spectra are 2204860, 1443843 and 499099, respectively 
(see Table 1). The spectra are classified as stars, galaxies, quasars and unknowns (mainly due to the poor 
quality of their spectra). DR1 includes both 1,487,200 spectra obtained during the first year of the general 
survey, and 717,660 spectra which were observed during the pilot survey (Luo et al. 2012). Due to the 
lack of network of photometric standard stars, the flux calibration is relative. The LAMOST spectra are 
unambiguously identified by a plate identifier, a spectrograph identifier and a fiber identifier (Song et al. 
2012; Ren et al. 2014). 

Based on the LAMOST data, lots of works have been done. For example, Ren et al. (2014) iden¬ 
tified 119 white dwarf-main sequence binaries, Zheng et al. (2014) reported a hypervelocity star with 
a heliocentric radial velocity of about 620 km/s, and Zhong et al. (2014) discovered 28 candidate 
high-velocity stars at heliocentric distances of less than 3 kpc. For more details see the website 
http://www.lamost.org/public/publication?locale=en. 

2.2 Optical spectra of symbiotic stars 

Optical spectra of known symbiotic stars have been catalogued by Allen (1984b) and Munari & Zwitter 
(2002). These spectra are characterized by the presence of absorption features typical of late-type giant 
stars (spectral type K or M). These include, among others. Cal, Fel, H 2 O, CO and TiO lines. Strong nebular 
Balmer line emission, as well as Hel, Hell and the [OIII], [Neill], [NV] and [FeVII] forbidden lines can 
also be seen. Additionally, more than 50% of symbiotic stars show broad emission features at A6830 A 
and a factor of ~ 4 weaker but similar emission feature at A7088 A, which only exists in symbiotic stars. 
The A6830 A band profile has a typical width of about 20 A and it often ranks among the 10 most intense 
lines in the optical region of symbiotic stars (which may reach 5% of the intensity of Ha). The AA6830 and 
AA7088 lines are due to Ramma scattering of the OVI resonance doublet AA1032, AA1038 by the neutral 
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V336 Car 



wavelength(A) 


Fig-1 V2416 Sgr (top panel) and V336 Car (bottom panel), an S-type and a D-type symbiotic 
star, respectively (Allen 1984a). The black and red lines are the spectra (continuum plus lines) 
and continua of the objects, respectively. We have normalized the flux into the range of [0, 1] 
following Eq. (2). 


In S-type stars (e.g. V2416 Sgr in Figure 1, top panel), the continuum has obvious absorption features 
originating from the red giant. Whilst the Balmer, Hel and Hell lines are strong, the forbidden lines such as 
[OIII], [Neill] are weak and even fade away. The Balmer decrement is generally steeper than Ho-/HA ~ 2.8 
(case B) (Kogure & Leung 2007). 

In D-type stars (e.g. V336 Car in Figure 1, bottom panel), the Balmer lines are strong and the forbidden 
lines are well developed in a wide range of excitation from [OI] up to [FeVII] and [NeV]. The Balmer 
decrement is also steeper than in the nebular case. The wide range of ionization levels in D-types reflects 
the fact that their nebulae are ionization-bounded whereas in at least some S-types the nebulae are density- 
bounded. 

D'-type stars (see Figure 6; note that this is one of our 4 identified LAMOST symbiotic binaries) have 
similar emission lines as for D-types. The most conspicuous property of D'-type stars is the presence of a 
G- or F-type giant instead of a cooler star (Miirset & Schmid 1999). 

According to above described features, we use the following criteria to classify an object as a symbiotic 
star (Mikolajewska et al. 1997; Mikolajewska et al. 2014; Belczynski et al. 2000): 

(a) Presence of the absorption features of a late-type giant: TiO, H 2 O, CO, CN and VO absorption 
bands, as well as Cal, Call, Fel and Nal lines. 

(b) Presence of strong HI and Hel emission lines and presence of emission lines from ions such as Hell, 
[OIII], [NV] and [FeVII] with an equivalent width exceeding 1 A. 
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LAMOST J220655.22-011747.5 



Fig. 2 LAMOST J220655.22-011747.5, is an MO type spectrum. The black line stands for the 
spectra of the object. The red line is the continuum fitted by the Global Continuum Fit method 
of Lee et al. (2008). 


2.3 Selection Method 


Based on the spectral properties of symbiotic stars and the feature of LAMOST data, we devise a selection 
method to find symbiotic stars. This follows a three-step procedure. 

(i) In order to directly compare the spectra, we first normalize the flux into the range of [0,1] following 


F n 


F\ — F min 

-Fmax -Fmin 


(1) 


where F max , F ru]u are the maximum and minimum fluxes of each spectrum, respectively. 

(ii) Unlike the vast majority of spectra observed by LAMOST, symbiotic stars display Balmer emission 
lines, especially the Ha and the H/3 lines. We therefore select symbiotic stars by analyzing the intensity of 
the Ha and H/3 lines, i.e. /(Ha) and /(H/3). These intensities are obtained from the following equation. 


n Ao+AA 

/(Ao) = / (-Fn - F c ) dX (2) 

where F c is the continuum fitted by the Global Continuum Fit method of Lee et al. (2008) (the 
[OIIIJA4340, A4959 and A5007 lines are removed when we fit the continuum). The AA for Ha, H/3 are 
set to 50 and 10 A, respectively. As the fluxes have been normalized, the intensity of the lines do not have 
truly physical significance. 

It is important to emphasize that fitting the continuum of M-type or later type spectra is subject to 
systematic uncertainties at the A5500 — 7500 A range, as shown in Figure 2. In order to efficiently select 
symbiotic stars, it becomes important to set proper values for the intensity of the Ha and H/3 lines. In our 
work, they were set to /(Ha) > 0.5 and /(H/3) > 0.2, which efficiently exclude spectra displaying the 
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Table 2 Spectra selected by our routine as symbiotic binaries 


Data 

Total 

Star 

M-type 

Step(ii) 

Step(iii) 

Symbiotic star 

DR1 

2,204,860 

1,944,406 

67,282 

44,983 

4,758 

1 

DR2 

1,443,843 

1,240,792 

54,888 

64,691 

7,631 

2 

DR3 

499,099 

437,299 

18,054 

15,408 

2,377 

1 


Notes: (1) Total stands for the number of the spectra we have used in each data release of the LAMOST. (2) 
Star, M-type columns stand for stars and M type spectra, respectively. (3) Step(ii) and Step(iii) columns stand for 
the selected out spectra by step(ii), step(iii) in Section 2.3. 


Table 3 Symbiotic stars observed by LAMOST 


Desig 

(LAMOST J) 

Date 

MJD 

IR type 

u 

9 

r 

i 

z 

Cool spectra 

DR 

Name 

122804.90-014825.7 

56592 

S 

15.88 

15.46 

12.18 

11.67 

11.92 

K3 

DR2 

— 

194957.58+461520.5 

56592 

S 

12.54 

— 

— 

— 

— 

M2 

DR2 

StHA 169 

202629.80+423652.0 

56946 

D 

15.05 

14.07 

13.49 

— 

— 

Mira 

DR3 

— 

214144.88+024354.4 

56204 

D' 

10.94 

10.33 

10.01 

11.39 

10.56 

G5 

DR1 

StHA 190 


Notes: Basic properties of the symbiotic stars, which are ordered by right ascension (R. A.) of the epoch J2000.0. 
The first and third symbiotic stars have been indentified by Downes & Keyes (1988), the others are new symbiotic 
stars. 


(iii) The Balmer series and other emission lines in symbiotic stars are very prominent, hence the contin¬ 
uum is generally flat in the optical spectrum (see Figure 1). We therefore average the normalized continuum 
F n in the 4000 — 7300 A wavelength range and exclude all spectra with the average value being above 0.4. 

After performing the above steps over the entire LAMOST spectroscopic database, we select 4758,7631 
and 2377 spectra in DR1, DR2, and DR3, respectively, as symbiotic star candidates (see Table 2). Visual 
inspection of the spectra revealed that most of these were planetary nebula and B[e] stars. Finally, we find 
four symbiotic stars, two of which have been catalogued in the past, the other two are new identifications 
(see Table 3). 


3 SYMBIOTIC STARS OBSERVED BY LAMOST 

Figures 3-6 show the LAMOST spectra of the 4 symbiotic stars we have identified. 

3.1 LAMOST J122804.90-014825.7 

The presence of strong Hell A4686 emission implies the presence of a hot component in the object 
(Figure 3). The spectrum is with a radial velocity of v r ~ (374 ± 2) km/s. And we can see the faint [FeVII] 
emission lines with the same Doppler shift in the inset of Figure 3. The spectrum displays also MgA5176 
and Call triplet absorption lines. It is classified as K3 type by the LAMSOT pipeline, which stands for the 
cool giant component of a symbiotic star. So we propose it is an S-type system. 

Based on the g, r and i magnitudes, the V magnitude is my ~ 13.5 mag 

(http://classic.sdss.org/dr4/algorithms/sdssUBVRITransform.html), and we assume a typical absolute 
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Fig. 3 LAMOST J122804.90-014825.7 is a new S-type symbiotic star, with a radial velocity 
i; r ~ (374 ± 2) km/s. The vertical dashed line in the inset stands for the wavelength of Hel 
A5876. The faint [FeVII] lines have the same velocity shift. F\ is relative flux. 


be about 5 kpc. The Galactic coordinates are l = 291.0062, b = 60.5226, so the vertical distance between 
the object and the Galactic disk is about 4.3 kpc. It seems to be a halo symbiotic star. 

3.2 LAMOST J194957.58+461520.5 

The LAMOST spectrum of J194957.58+461520.5 is very similar to the one presented by Downes & Keyes 
(1988) for the 4000 - 7300 A range. LAMOST J194957.58+461520.5 is included in the Belczynski et al. 
(2000) catalogue of symbiotic stars as StHA 169. The spectrum shows a strong and broad Ha emission line. 
The other Balmer lines, as well as Hel and HeIIA4686 lines are also clearly present. At the same time there 
is a faint feature at AA6830 as shown in Figure 4. The spectrum displays obvious TiO absorption bands 
and Call triplet absorption lines, and the star is classified as M2 by the LAMOST pipeline. It is an S-type 
symbiotic star. 

3.3 LAMOST J20262980+4236520 

The spectrum of this object is shown in Figure 5. The Balmer and [OIII] lines are clearly present. Weak 
[Nil] A6548, A6584 lines at the both sides of the Ha line can also be seen, however, the Hel and Hell lines 
are very weak. There is an obvious AA6830 line (/ (AA6830) // (Ha) ~ 7.6%) and a AA7088 line. The flat 
continuum, and the intensity ratio I ([OIII]A5007) /I (H/3) ~ 0.89 and I ([OIII]A4363) // (H 7 ) ~ 4.29 
(Gutierrez-Moreno & Moreno 1995; Mikolajewska et al. 2014; Rodrfguez-Flores et al. 2014) suggests that 
LAMOST J20262980+4236520is a D-type symbiotic star, so its cool giant is a variable. The signal to noise 
ratio (S/N) is relatively low at the iz bands where the spectrum is also affected by systematic uncertainties 
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165 LAMOSTJ 194957.58+461520.5 



Fig. 4 LAMOST J1194957.58+461520.5 is a S-type symbiotic star, which has been observed 
by Downes & Keyes (1988), and listed as StHA 169 in Belczynski et al. (2000) catalogue of 
symbiotic stars. F\ is relative flux. 


LAMOST J202629.80+423652.0 



Fig. 5 LAMOST J20262980+4236520 is a new D-type symbiotic star. The red part of the spec¬ 
trum is affected by systematic uncertainties due to sky subtraction. F\ is relative flux. 

3.4 LAMOST J214144.88+024354.4 

LAMOST J214144.88+024354.4 is a D'-type symbiotic star. It is included in the catalogue by 
Belczynski et al. (2000) as StHA190. The spectrum, shown in Figure 6, is nearly identical to the one pre¬ 
sented by Downes & Keyes (1988) in the 4000 — 7300 A wavelength range. The spectrum is classified as 
G5 type by the LAMOST pipeline, which is consistent with the D'-type nature of this symbiotic binary. 
High-resolution optical spectra of StHA 190 were analyzed by Smith et al. (2001), who derived the giant 
temperature, T e g = 5300 ± 150 K and [Fe/H] ~ 0 based on fitting the spectrum with model atmospheres. 
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Fig. 6 LAMOST J214144.88+024354.4 is a D'-type symbiotic star, included in the 
Belczynski et al. (2000) list as StHA 190. F\ is relative flux. 

4 DISCUSSION 

We used 4,147,802 LAMOST spectra to find symbiotic stars, of which 3,622,479 are classified as stars by 
the LAMOST pipeline (see Table 2). Since about 1% of the stars are red giants in our solar neighborhood 
(http://pages.uoregon.edu/jimbrau/astrl22/Notes/Chapterl7.html), a rough estimate of the number of the 
red giant stars in these spectra is around 30,000. Magrini et al. (2003) proposed that the number of symbiotic 
stars can be taken to be 0.5% of the total number of red giants. Therefore there should be more than 100 
symbiotic stars to be selected out. But we should note that the estimated number of Galactic symbiotic stars 
from Magrini et al. (2003) is 1 -2 orders of magnitude higher than that given by others as shown in Section 1. 
So the real symbiotic stars among the LAMOST stellar spectra may be less than 100. As the first step, we 
only identified 4 symbiotic binaries, far below what we expected. By comparing the coordinates between 
the known symbiotic stars listed in the catalogue of Belczynski et al. (2000) and the footprint of LAMOST 
DR1, there should be 9 known symbiotic stars in the area of the footprint (see Table 4 ). According to 
the plate design of the LAMOST, for brighter stars (r < 14 mag), only those with 9 < J < 12.5 mag 
from the 2MASS catalogs are chosen as potential targets (Luo et al. 2015). It immediately rules out 8 of 
the 9 symbiotic stars. For the fainter stars (r > 14 mag), there are 5 conditions to select potential targets 
(Yuan et al. 2015), which might exclude the symbiotic star Draco C-l. The StHA 190 has not been excluded, 
probably because its J-band magnitude is very close to the target selection condition of brighter stars. In the 
future work, we will apply the selection method to the SDSS (the Sloan Digital Sky Survey) released data 
to find more symbiotic stars and check its efficiency. 

5 CONCLUSION 

We have found 4 symbiotic stars among 4,147,802 spectra released by the LAMOST DR1, DR2 and DR3. 
The objects LAMOST J194957.58+461520.5 and LAMOST J214144.88+024354.4 are included in the 
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Table 4 Symbiotic stars in the footprint of LAMOST DR1 



Name 

RA 

Dec 

J 

r 


EG And 

00 44 37.1 

+40 40 45.7 

3.65 

6.6 

BD Cam 

03 42 09.3 

+63 13 00.15 

1.31 

4.5 

UV Aur 

05 21 48.8 

+32 30 43.1 

4.03 

8.3 

TX CVn 

12 44 42.0 

+36 45 50.6 

7.47 

8.7 

TCrB 

15 59 30.1 

+25 55 12.6 

5.70 

2.1 

AGDra 

16 01 40.5 

+66 48 09.5 

7.15 

9.2 

Draco C-1 

17 19 57.6 

+57 50 04.9 

14.38 

16.7 

StHA 190 

21 41 44.8 

+02 43 54.4 

8.71 

10.0 

AG Peg 

21 51 01.9 

+ 12 37 29.4 

5.00 

5.8 


Notes: The symbiotic stars are in the area of the footprint of LAMOST DR1. RA and Dec are the celestial 
coordinates. 


can be selected out by our method. The LAMOST spectra of these two objects are nearly identical to those 
provided by Downes & Keyes (1988). LAMOST J12280490-014825.7 is a new symbiotic star and is con¬ 
sidered to be a Halo S-type with a radial velocity of v r ~ (374 ± 2) km/s. LAMOST J20262980+4236520 
is also a new symbiotic star classified as D-type. With the operating of the LAMOST, more and more 
symbiotic star spectra will be obtained. 
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